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Abstract. The escape probability is a deterministic concept that quantifies some as- 
pects of stochastic dynamics. This issue has been investigated previously for dynami- 
cal systems driven by Gaussian Brownian motions. The present work considers escape 
probabilities for dynamical systems driven by non-Gaussian Levy motions, especially 
symmetric a-stable Levy motions. The escape probabilities are characterized as solu- 
tions of the Balayagc-Dirichlct problems of certain partial differential-integral equations. 
Differences between escape probabilities for dynamical systems driven by Gaussian and 
non-Gaussian noises arc highlighted. In certain special cases, analytic results for escape 
probabilities are given. 
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1. Introduction 

Stochastic dynamical systems arise as mathematical models for complex phenomena 
in biological, geophysical, physical and chemical sciences, under random fluctuations. A 
specific orbit (or trajectory) for such a system could vary wildly from one realization to 
another, unlike the situation for deterministic dynamical systems. It is desirable to have 
different concepts for quantifying stochastic dynamical behaviors. The escape probability 
is such a concept. 

Brownian motions are Gaussian stochastic processes and thus are appropriate for mod- 
eling Gaussian random fluctuations. Almost all sample paths of Brownian motions are 
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continuous in time. For a dynamical system driven by Brownian motions, almost all or- 
bits (or paths or trajectories) are thus continuous in time. The escape probability is the 
likelihood that an orbit, starting inside an open domain D, exits this domain first through 
a specific part T of the boundary dD. This concept helps understand various phenomena 
in sciences. One example is in molecular genetics [2B]. The frequency of collisions of two 
single strands of long helical DNA molecules that leads to a double-stranded molecule is 
of interest and can be computed by virtu of solving an escape probability problem. It 
turns out that the escape probability satisfies an elliptic partial differential equation with 
properly chosen boundary conditions [261 CHI [25j IS]- 

Non-Gaussian random fluctuations are widely observed in various areas such as physics, 
biology, seismology, electrical engineering and finance [29| IT6| 120] . Levy motions are a 
large class of non-Gaussian stochastic processes whose sample paths are discontinuous 
in time. For a dynamical system driven by Levy motions, almost all the orbits X t are 
discontinuous in time. In fact, these orbits are cadlag (right-continuous with left limit at 
each time instant), i.e., each of these orbits has countable jumps in time. Due to these 
jumps, an orbit could escape an open domain without passing through its boundary. In 
this case, the escape probability is the likelihood that an orbit, starting inside an open 
domain D, exits this domain first by landing in a target domain U in D c (the complement 
of domain D). 

As we see, the escape probability is defined slightly differently for dynamical systems 
driven by Gaussian or non-Gaussian processes. Although the escape probability for the 
former has been investigated extensively, the characterization for the escape probability for 
the latter has not been well documented as a dynamical systems analysis tool for applied 
mathematics and science communities. See our recent works [HI EH] for numerical analysis 
of escape probability and mean exit time for dynamical systems driven by symmetric 
a-stable Levy motions. 

In this paper, we carefully derive a partial differential-integral equation to be satisfied 
by the escape probability for a class of dynamical systems driven by Levy motions, es- 
pecially symmetric a-stable Levy motions. Namely the escape probability is a solution 
of a nonlocal differential equation. We highlight the differences between escape probabil- 
ities for dynamical systems driven by Gaussian and non-Gaussian processes. These are 
illustrated in a few examples. 

More precisely, let {X t , t ^ 0} be a Revalued Markov process defined on a complete 
filtered probability space (f2, J 7 , {J-"t}t^o, P)- Let D be an open domain in M. d . Define the 
exit time 

T D c := inf{t > : X t E D c }, 

where D c is the complement of D in M. d . Namely, is the first time when X t hits D c . 

When X t has almost surely continuous paths, i.e., X t is either a Brownian motion or a 
solution process for a dynamical system driven by Brownian motions, a path starting at 
x G D will hit D° by hitting dD first (assume for the moment that dD is smooth). Thus 
t~d c = T dD- Let T be a subset of the boundary dD. The likelihood that X t , starting at 
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x, exits from D first through T is called the escape probability from D to T, denoted as 
p(x). That is, 

P (x) = ¥{x TaD e r}. 

We will verify that (Section 13. 2p the escape probability p(x) solves the following Dirichlet 
boundary value problem: 

Cp = 0, x E D, ,^ 
PldD = *P, 

where £ is the infinitesimal generator of the process X t and the boundary data ip is defined 
as follows 



ip(x) 



l, xeT, 

0, xedD\T. 



When X t has cadlag paths which have countable jumps in time, i.e., X t could be either 
a Levy motion or a solution process of a dynamical system driven by Levy motions, the 
first hitting of D c may occur somewhere in D c . For this reason, we take a subset U of 
D c , and define the likelihood that X t exits firstly from D by landing in the target set U 
as the escape probability from D to U, also denoted by p(x). That is, 

p(x) = F{X TDC G U}. 

We will demonstrate that (Section 13 .4j) the escape probability p(x) solves the following 
Balayage-Dirichlet boundary value problem: 

Ap = 0, xeD, 
p\d» = 

where A is the characteristic operator of X t and ip is defined as follows 

( v _ / l, xeU, 
kp[x> ~ \ o, x e D c \ u. 

Therefore by solving a deterministic boundary value problem (pQ) or ([2]), we obtain the 
escape probability p(x). 

This paper is arranged as follows. In Section^ we introduce Balayage-Dirichlet problem 
for discontinuous Markov processes, and also define Levy motions. The main result is 
stated and proved in Section [31 In Section HI we present analytic solutions for escape 
probabilities in a few special cases. 

2. Preliminaries 

In this section, we recall basic concepts and results that will be needed throughout the 
paper. 



2.1. Balayage-Dirichlet problem for discontinuous Markov processes. The fol- 
lowing materials are from [3J [T7J |27J [7J [191 H2]- Let Q be a locally compact space with a 
countable base and £f be the Borel cx-field of Also, ^ is adjoined to Q as the point at 
infinity if Q is noncompact, and as an isolated point if Q is compact. Furthermore, let ^ 
be the cx-field of Borel sets of Q q = Q U {c}. 

Definition 2.1. A Markov process Y with state space (G,&) is called a Hunt process 
provided: 

(i) The paths functions t — > Y t are right continuous on [0, oo) and have left-hand limits 
on [0, () almost surely, where ( := inf{t : Y t = c}. 

(ii) Y is strong Markov. 

(Hi) Y is quasi-left- continuous: whenever {r n } is an increasing sequence of Ft- stopping 
times with limit r , then almost surely Y Tn — > Y T on {r < oo}. 

Definition 2.2. Let G be an open subset of Q and Y t (x) be a Hunt process starting at 
x G G. A nonnegative function h defined on Q is said to be harmonic with respect to Y t 
in G if for every compact set K C G, 

E[h(Y TKC (x))\ = h{x), xeG. 

Definition 2.3. Let f be nonnegative on G c . We say h defined on Q solves the Balayage- 
Dirichlet problem for G with "boundary value" f , denoted by (G, f), if h = f on G c , h is 
harmonic with respect to Y t in G and further satisfies the following boundary condition: 

Wz e dG, h(y) — > f(z), as y — > z from inside G. 

A point z G dG is called regular for G c with respect to Y t (z) if 

P{r G c = 0} = 1. 

Here G is said to be regular if any z G dG is regular for G c . 

Let p be a metric on Q compatible with the given topology. Let Xq be the family of 
functions g ^ bounded on Q and lower semicontinuous in G such that Va; G G, there is 
a number Ag{x) satisfying 

Eb( v;,W)]- gM ^ 

E[t £ \ 

where r e := inf{t > : p(Y t (x),x) > e}. We call A with domain Xq the characteristic 
operator of Y t relative to G. If L with domain T>q is the infinitesimal generator of Y t 
relative to G, V G C X G and 

Af = Cf, feV G . 

(c.f. nm) 

We quote the following result about the existence and regularity of the solution for the 
Balayage-Dirichlet problem. 

Theorem 2.4. (^7\) 

Suppose that G is relatively compact and regular, and f is nonnegative and bounded on 
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G c . If f is continuous at any z G dG, then h(x) = K[f(Y TGC (x))} is the unique solution 
to the Balayage-Dirichlet problem (G, f), and Ah(x) = for h G Ig- 

2.2. Levy motions. 

Definition 2.5. A process Lt, with Lq — a.s. is a d- dimensional Levy process or Levy 
motion if 

(i) L t has independent increments; that is, L t — L s is independent of L v — L u if (u, v) D 
(M) = 0; 

(ii) L t has stationary increments; that is, L t — L s has the same distribution as L v — L u 
if t — s = v — u > 0; 

(Hi) L t is stochastically continuous; 
(iv) L t is right continuous with left limit. 

The characteristic function for L t is given by 

E (exp{i(z, L t )}) = exp{t^{z)}, z G M d , 

where (•, •) is the scalar product in M d . The function \I/ : M. d — > C is called the charac- 
teristic exponent of the Levy process L t . By the Levy-Khintchine formula, there exist a 
nonnegative-definite d x d matrix Q, a measure v on IR d satisfying 

z/({0}) = 0and / (|u| 2 A l)v(du) < oo, 

JR d \{0} 

and 7 G M d such that 

V{z) = i(z,7> ~-(z,Qz) 

+ [ (e^ u >-l-i(z, U )l M<1 )Kdw). (3) 

*/R d \{0} 

The measure v is called the Levy measure of L t , Q is the diffusion matrix, and 7 is the 
drift vector. 

We now introduce a special class of Levy motions, i.e., the symmetric a-stable Levy 
motions Lf . 

Definition 2.6. For a G (0, 2) . A d-dimensional symmetric a-stable Levy motion Lf is 
a Levy process with characteristic exponent 

V{z) = -C\z\ a , z G R d , (4) 

where 

_ 1/2 r((l + q)/2)r(d/2) 
r((d+a)/2) ' 

(c.f. pU Page 115] for the above formula of C .) 

5 



Thus, for a d-dimensional symmetric a-stable Levy motion Lf, the diffusion matrix 
Q = 0, the drift vector 7 = 0, and the Levy measure v is given by 

u(du) = . du, 

V I \ u \d+<* 

where 

aT((d + a)/2) 



2i-«7r d / 2 r(l -a/2)' 

(c.f. Page 1312] for the above formula of C^ a -) Moreover, comparing (j3J) with ([3]), we 
obtain 

—C\z\ a = I {e^ -l-i{z,u)l H ^)^du. 

jR d \{0} \ u \ 

Let Co(M. d ) be the space of continuous functions / on M d satisfying lim f(x) = with 

\x\— >oo 

norm ||/||c (R d ) = SU P 1/(^)1- Let C$(M. d ) be the set of / G C (M d ) such that / is two times 



different iable and the first and second order partial derivatives of / belong to Co(M d ). Let 
C^°(R d ) stand for the space of all infinitely differentiate functions on M. d with compact 
supports. Define 

(C a f)(x) : = f (f(x + u)- f(x) - (d x f(x), m)1|«|<i) ■ \Xa du 



on C$(R d ). And then for (6R d 

(C a e^)(x) = e«**> [ (e<™ - 1 - i(^>l H<1 )-%^. 

•/R d \{0} \ u \ 

By Courrege's second theorem ([H Theorem 3.5.5, p. 183]), for every / G C^°( 
(£«/)(*) = j^-J e^[e-^ z \C a e^){x)}f{z)dz 



co Cud\ 



1 



,i{z,x) 

^\z\ a f{z)dz 



(e^-l-il^)!^)^ 



f(z)dz 



(2vr) d / 2 J Rd [jRd\ {0} \u\ 

C 

~ (2n) d / 2 

= c-[-{-AY' 2 f]{x). 

Set p t := L t — L t -. Then p t defines a stationary (J^-adapted Poisson point process 
with values in M. d \ {0} ([12]). And the characteristic measure of p is the Levy measure 
v. Let N p ((0, t], du) be the counting measure of p t , i.e., for B G B(M d \ {0}) 

N p ((Q,t],B):=#{0<s^t:p s eB}, 

where # denotes the cardinality of a set. The compensator measure of N p is given by 

N p {(0,t],du) := N p ({0,i\,du)-tv(du). 



The Levy-Ito theorem states that for a symmetric ct-stable process L t , 

(i) for 1 ^ a < 2, 

L t = f f uN p (ds,du)+ f f uN p (ds,du), 
Jo J\u\<:i Jo J\u\>l 

(ii) for < a < 1, 

L t = / / uN p (ds,du). 
Jo Jm d \{o} 

3. Boundary value problems for escape probability 

In this section, we formulate boundary value problems for the escape probability as- 
sociated with Brownian motions, SDEs driven by Brownian motions, Levy motions and 
SDEs driven by Levy motions. For Levy motions, in particular, we consider symmetric 
a-stable Levy motions. We will see that the escape probability can be found by solving 
deterministic partial differential equations or partial differential-integral equations, with 
properly chosen boundary conditions. 



J I L 



D 



x — 5 x x + 5 



Figure 1 . A particle executing unbiased random walk in a bounded interval 

3.1. Boundary value problem for escape probability of Brownian motions. Sup- 
pose that a particle executes an unbiased random walk on a straight line. Let D = (a, b). 
Figure [1] shows the random walk scenario. That is, a particle moves according to the 
following rules (|18j): 

(i) During the passage of a certain fixed time interval, a particle takes 1 step of a certain 
fixed length § along the x axis. 

(ii) It is equally probable that the step is to the right or to the left. 

If the particle starting from x G D eventually escapes D by crossing the boundary b, 
then it must have moved to one of the two points adjacent to x first and then crossed the 
boundary. Thus 

P( x ) = ^[p(x - 5) + p(x + 5)], 
for x G D. By Taylor expansion on the right hand side to the second order, we have 

l -p'\x) = 0. 

The boundary conditions are 

limp(x) = 1, \imp(x) = 0, 

x—>b x— >a 
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since the nearer the particle starts to b, the more likely it will first cross the boundary 
through b. 

Note that the limit of the random walk is a standard Brownian motion Wt, that is, 

(i) W has independent increments; 

(ii) for < s < t, Wt — W s is a Gaussian random variable with mean zero and variance 
(t-s). 

Thus, the escape probability p(x) of a standard Brownian motion from D through the 
boundary b satisfies 

f Ap(z) = 0, 
p(b) = 1, 
p(a) = 0, 

where ~A = \d xx is the infinitesimal generator for a scalar standard Brownian motion 
W t . 

3.2. Boundary value problem for escape probability of SDEs driven by Brow- 
nian motions. Some results in this subsection can be found in [2TJ Chapter 9]. 

Let {W(t)} t ^o be an m-dimensional standard J-"radapted Brownian motion. Consider 
the following stochastic differential equation (SDE) in M. d : 

X t (x) = x+ [ b(X s (x))ds+ [ a(X s (x))dW s . (5) 
Jo Jo 

We make the following assumptions about the drift b : M. d i— y M. d and the diffusion coeffi- 
cient a : R d R d x R m . 

\b(x)-b(y)\^X(\x-y\), 
\a(x) - a(y)\ ^ 7(|x -y\). 

Here A and 7 are increasing concave functions with the properties A(0) = 7(0) = 0, 
fn^ Trrdw = L . -^U-dw = 00. 

J0+ X(u) J0+ y (u) 

Under (H£ a ), it is well known that there exists a unique strong solution to Eq.fl5|)([30]). 
This solution is denoted by X t (x). 

We also make the following assumption. 

(H^.) There exists a £ > such that for any x,y e D 

(y,aa*{x)y) ^ £\y\ 2 . 

This condition guarantees that the infinitesimal generator 

^ ^ 0^ 

i=l 1 i,j=l 1 J 

for Eq.(E]) is uniformly elliptic in D, since then the eigenvalues of 00* are away from in 
D. Here the matrix [a^] := ^a(x)a*(x). 



Figure 2. Escape probability for SDEs driven by Brownian motions: An 
annular open domain D with a subset T of its boundary dD 



Let D be an open annular domain as in Figure |2j In one dimensional case, it is just an 
open interval. Let T be its inner (or outer) boundary. Taking 



This means that, for this specific if;, E,[ijj(X TgD (x))} is the escape probability p(x), which 
we are looking for. 

We need to use [2TJ Theorem 9.2.14] or [9] in order to see that the escape probability 
p(x) is closely related to a harmonic function with respect to X t . This requires that the 
boundary data ip to be bounded and continuous on dD. For the domain D taken as 
in Figure [21 with T the inner boundary (or outer) boundary, the above chosen ip in ([6]) 
is indeed bounded and continuous on 3D. Thus, we have the following result by [2Tj 
Theorem 9.2.14]. 

Theorem 3.1. The escape probability p(x) from an open annular domain D to its inner 
(or outer) boundary Y, for the dynamical system driven by Brownian motions is the 
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x e r, 
ie3D\r, 



(6) 



we have 




F{u:X TdD (x)eT} 
p(x). 



solution to the following Dirichlet boundary value problem 

Cp = 0, 
P\r = 1, 
p\dD\r = 0. 




Figure 3. A particle executing Levy motion in a bounded interval 

3.3. Boundary value problem for escape probability of symmetric a-stable Levy 
motions. Assume a particle is taking a one- dimensional Levy flight, where the distribu- 
tion of step sizes is a symmetric a-stable distribution (Figure |3]). Let p(x) denote the 
escape probability of the particle starting at x in D = (a, b) and then first escapes D over 
the right boundary b. It could first move to somewhere inside D, say x + y G D, and then 
achieve its goal by jumping over the right boundary b from the new starting point x + y. 
More precisely, 

p(x) = / P{the first step length is y}p(x + y)dy. (7) 

jR\|0j 



According to [2], the symmetric a-stable probability density function is the following 



faM 



-^i^^sin^f -«argy)], < a < 1, 



^Er=o(-l) fe£ ^ fc cos[Mi)L 1 < a <2, 

where argy = ir when y < 0. 

For < a < 2, the asymptotic expansion has also been given by [2] as follows 



■i ELi ^'W 1 sin ~ a ^v)] + odyl"* 1 ^- 1 ), Is/1 °°, 
iELo(-i)*n#Vcos [*(=)] + (|i/r 1 ), |y| o, 

'CxCaVlj/l^ + odj/l- 2 "- 1 ), M-)-oo, 
C 2 (a) + o(|y| 2 ), 0, 



where Ci(a) = - sin(^)r(l + a) and C 2 (a) = - r ^ a ^ . Take A > large enough and fix 
it. Thus, 



= / f afi (y)\p(x + y) - p(x)]dy 

fa,o(y)[p( x + y) -p( x )]dy 



(-N,N)\{0} 
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+ 

Ir\(-n,n) 
=: h + h- 

For Ii, by self-affine property in [28], we obtain 

II 



h 



f a ,o(N ■ —)[p(x + y) - p(x)]dy 

{-N,N)\{0} iV 

Lo,w^ b(i+ ^ p(i)ld ' 



C x {a) 



(-n,n U o } N^(\y\/Ny+« 
f Ci(a 

(-N,N)\{0} 



\p(x + y) -p(x)]dy 



\ y \l+a 



\p(x + y) -p(x))dy. 



For I 2 , we calculate 



-A' 



f a ,o(y)[L - P(x)]dy + / f a ,o(y)[0 -p(x)]dy 



00 



A 



Ci(g) 



-A" 



;1+2q 



Ci(g) 



[1 -p(ar)]dj/ 
1 



-2a 



p(x)d?/ 



00 



jv y 



l+a 



dy 



CM) 



M\[-iV,Ar] 



Ci(a)p(x + y) 



I l+a 



dy 



(9) 



Jr\[-n,n] \y\^" Jr\[-n,n] 

Note that for < a < 1, by the fact that the integral of an odd function on a symmetric 
interval is zero, it holds that 



/ p\ x )yh\ymrAi~k dy = p\ x ) c M) [ 
Jr\{o\ \y\ J{\i 



v 



Aiwi<i}\{o} \y\ 

Thus, putting (jSJ), (jSJ) and ({TO]) together, we have for < a < 1 

/ [p(x + y)- p(x) - p'(x)yl { \ y]<1} ] ^^dy = 0. 

Jm\{o} \y\ 

Moreover, C\(a) = C\ )OL . 

For a G [1, 2), we only divide l\ into two parts l\\ and /12, where 



Ti-dy = 0. 

l+a ^ 



(10) 



In : = 



{|vl<e}\{0} 



/ a ,o(y)[p(z + y) -p(^)]dy, 
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fa,o{y)\p{x + y)-p(x)]dy, 

(-N,N)\{-q,q) 

and q > is a small enough constant. 

For In, by Taylor expansion and self-affine property in [26], we get 

fa,o(y)\p(x + y) -p(x)]dy 

{\y\<s}\{o} 

f a ,o(- ■ Qy)p\ x )y^y 

{|y|<e}\{o} Q 

fafi(g) u , , 

V {x)ydy 



{\y\<?}\{0} (e\y\) 1+a 

CAa 



{ivi«>\{o> (^) 1+Q • (<?lvl) 1+a 

/ v ( x )yh\y\^}rji^ d y 
jR\io\ \y\ 



p'{x)ydy 



For I12, we apply the same technique as that in dealing with Ii for a G (0, 1). 
Next, by the similar calculation to that for a G (0, 1), we obtain for a G [1,2) 

[p(x + y)- p(x) - p'(x)yI M< i } ] C ^l^ V = 0. 

m\{o} \y\ 

Since the limit of the Levy flight is a symmetric a-stable Levy motion L", the escape 
probability p(x) of a symmetric a-stable Levy motion, from D to [b, oo) satisfies 

-(-A)fp(x) = 0, 
p(x)\[b,oo) = 1, 

P(z)|(-QO,a] = 0. 

Note that — (— A)§ is the infinitesimal generator for a scalar symmetric a-stable Levy 
motion Lf. 

3.4. Boundary value problem for escape probability of SDEs driven by general 
Levy motions. Let L t be a Levy process independent of W t . Consider the following SDE 
in R d 

X t {x) = x+ [ b(X s {x))ds+ [ a(X s (x))dW s + L t . (11) 
Jo Jo 

Assume that the drift b and the diffusion a satisfy the following conditions: 

(Hb) there exists a constant Cf, > such that for x, y G W 1 

\b(x) - b(y)\ ^ C b \x - y\ ■ logflx - y]' 1 + e); 

(H CT ) there exists a constant C a > such that for x, y G W 1 

\a(x) - a(y)\ 2 < C a \x - y\ 2 ■ log(\x - y^ 1 + e). 
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Under (H&) and (H CT ), it is well known that there exists a unique strong solution to 
Eq. ( ITT]) (see [23]). This solution will be denoted by X t (x). Moreover, X t (x) is continuous 
in x. 

Lemma 3.2. The solution process X t (x) of the SDE ( ITT]) is a strong Markov process. 

Proof. Let 77 be a (J-t)^o-stopping time. Set 

Qt := a{W v+t - Wr,, L v+t -L v }UAf, t> 0, 

where Af is of all P-zero sets. That is, Q t is a completed a-algebra generated by W v+t — W v 
and L v+t — L v . Besides, Q t is independent of T t - Let X(x,r},T] + t) denote the unique 
solution of the the following SDE 

rn+t rv+t 
X(x : r) : r) + t) =x+ b(X(x,r),s))ds+ a(X(x, 77, s)) dW s + L v+t — L v . (12) 

J rj J T) 

Moreover, X(x,r],r] + 1) is (^-measurable and X(x,0, t) = X t (x). By the uniqueness of 
the solution to (fT2]) . we have 

X(x, 0, 77 + t) = X(X(x, 0, 77), 77, 77 + 1), a.s.. 

For any bounded measurable function g, 

E^^+^x))!^] = E[^(X(x,0,77 + t))|^] 

= E[g(X(X(x,0,r l ), V , V + t))\J v } 

= ^[g( x (y,V:V + t))]\y=x(x,o, V ) 

= E[g(X(y,0,t))]\ y=x{xM . (13) 

Here the last equality holds because the distribution of X(y, 77, 77 + t) is the same to that 
of X(y,0,t). The proof is completed since (fTB"]) implies that 

E[g(X v+t (x))n = ElgiX^x^X.ix)]. 

□ 

Because L t has cadlag and quasi-left-continuous paths([21]), X t (x) also has cadlag 
and quasi-left-continuous paths. Thus by Lemma 13.21 and Definition 12.11 above, we see 
that X t (x) is a Hunt process. Let D be a relatively compact and regular open domain 
(Figure H] or Figure [5]). Theorem 12.41 implies that E[ip(X TDC (x))} is the unique solution 
to the Balayage-Dirichlet problem (D,ip), under the condition that ip is nonnegative and 
bounded on D c . Set 



ip(x) 



1, x e U, 
0, xeD c \U. 

Then ip is nonnegative and bounded on D c . We observe that 



E[cp(X TDC (x))} = / cp(X TDC (x))d¥(cu) 

<{W.X TDC {X)£U} 

13 



Figure 4. Escape probability for SDEs driven by Levy motions: an open 
annular domain D, with its inner part U (which is in D c ) as a target domain 




Figure 5. Escape probability for SDEs driven by Levy motions: a general 
open domain D, with a target domain U in D c 

+ / if(X TDC (x))d¥(cu) 
J {w.x TDC (x)eD-\u} 
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= p(x). 



This means that, for this specific cp, E[cp(X Tnc (x))] is the escape probability p{x) that we 
are looking for. By the definition of the characteristic operator, p6/fl and by Theorem 
12.41 Ap(x) = 0. Thus we obtain the following theorem. 

Theorem 3.3. Let D be a relatively compact and regular open domain, and let U be a set 
in D c . Then the escape probability p(x), for the dynamical system driven by Levy motions 
fTTTj) . from D to U , is the solution of the following Balayage-Dirichlet problem 

Ap = 0, 
P\u = 1, 
p\d°\u = 0> 

where A is the characteristic operator for this system. 

Remark 3.4. Unlike the SDEs driven by Brownian motions, a typical open domain D 
here could be a quite general open domain (Figure [3J), as well as an annular domain 
(Figure^. This is due to the jumping properties of the solution paths. It is also due 
to the fact that, in Theorem \2.4\ the function f is only required to be continuous on the 



boundary 3D (not on the domain D°). 



Finally we consider the representation of the characteristic operator A, for a SDE driven 
by a symmetric a-stable Levy process Lf , with a G (0, 2): 

X t (x)=x + [ b(X s (x))ds+ [ <r(X a (x))dW a + Lf. (14) 
Jo Jo 

Let us first consider the case of 1 ^ a < 2. For / G Cg(lR d ), applying the Ito formula 
to f(X Te (x)), we obtain 

f(X Ts (x)) - f(x) = f E (d y f(X s ),b(X s ))ds + r (d y f(X s ),a(X s )dW s ) 

Jo Jo 

+ r f (f(X s + u)- f(X s )) N p (ds, du) 

Jo J\u\^l 



P [ (f(X s + u)-f(X s ))N p (ds,du) 

Jo J\u\>l 

f(X s )) a ik (X s )a kj (X s )ds 



>0 J\u\>l 

i n ( d 2 



2 Jo \dyidyj 
■ P I (f(X. + u) - f(X s ) 

Jo J\u\^l 

C 

-{dyf{X s ),u))—-^duds. 
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Here and hereafter, we use the convention that repeated indices imply summation from 1 
to d. Taking expectation on both sides, we get 



E[f(X Te (x))]-f(x) 



E / (d y f(X s ),b(X s ))ds + -E 



d 2 



-E 



o \dyidyj 
(f(X s + u)-f(X s )-(d y f(X s ),u)) u 



f(X s ) ) a ik (X s )a kj (X s )ds 



duds. 



\u\ 



'0 JM d \{0} 

The infinitesimal generator £ of Eq. ffll]) is as follows p] 

1 / d 2 

(Cf)(x) = (d x f(x),b(x)) + - ( ^ Q „ f(x) ) a ik (x)a kj (x) 



2 \dxidxj 



(f( x + u )-f(x)-(d x f(x),u)) 



Cd,a 



d+a 



du. 



So, 



E[f(X Te (x))]-f(x) 



E 



E£(Cf)(X s )ds Ej^(Cf)(x)ds 



E\r P 



E\r F 



E£\(Cf)(X s )-(Cf)(x)\ds 



El 



7V 



^ sup \(Cf)(y) - (Cf)(x)\. 

\y-x\<£ 



Because is continuous in x, 

E[/(X Te (*))]-/(x) 



Af{x) = lim 

e|0 



Eh, 



Similarly, we also have A = C for < a < 1. 



Remark 3.5. The above deduction tells us Af = Cf for f G Co(IR d ). If the considered 
driving process is not a symmetric a-stable Levy motion, the domain of C is unclear and 
thus A = C may not be true. The corresponding escape probability p(x) is the solution of 
the following Balayage-Dirichlet problem (in terms of operator C, instead of A): 

Cp = 0, 
P\u = 1, 
pW\u = o. 



4. Examples 



In this section we consider a few examples. 
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Figure 6. Escape probability for 1-dimensional Brownian motion in Ex- 
ample 14.11 r = 2 

Example 4.1. In 1-dimensional case, take D = (— r, r) and F = {r}. For each x G D, 
the escape probability p(x) of X t = x + W t from D to T satisfies the following differential 
equation 

{\p"(x) = 0, x E (— r, r), 
p(r) = 1, 
p(—r) = 0. 

We obtain that p(x) = f or x £ [_ r? r j_ jf i s a straight line (See Figure\M). 

In 2-dimensional case, take D = {x G lR 2 ;r < \x\ < R} and T = {x G 1R 2 ; \x\ = r}. 
For every x G D, the escape probability p(x) of X t = x + Wt from D to T satisfies the 
following elliptic partial differential equation 

lAp(x) = 0, x G D, 

p{x)\\x\=r = 1, 
p(x)\\x\=r = 0. 

By solving this equation, we obtain that p(x) = It is plotted in Figure^ 

Example 4.2. Consider the following SDE driven by Brownian motions: 

dX t = b(X t )dt + a(X t )dW t , 

where b and (non-zero) a are real functions. When b and o satisfy (H^ CT ), the equation 
has a unique solution which is denoted as X t . We take D = (— r, r) and T = {r}. For 

17 



Figure 7. Escape probability for 2-dimensional Brownian motion in Ex- 
ample ED r = 2, R = 4 



each x G D, under the condition (H 2 ), the escape probability p(x) satisfies 

^a 2 (x)p" (x) + b(x)p (x) = 0, x G (— r, r), 



p(r) = 1, 
p(—r) = 0. 



TTie solution is 



\ e ( ) dy 

P x = TT\ 

J e J - r ^) dy 



for x G [— r, r] . See Figure 



Example 4.3. In 1- dimensional case, take D = (— r, r) and U = [r, oo). For eac/i x £ D 
and a symmetric a-stable Levy process L%, the escape probability p(x) of X t = x + 
from D to U satisfies the following differential-integral equation 

— (— A)fp(x) =0, x G (— r, r), 

p{x)\[r,oo) = 1, 
p(ar) I (_«,,_,.] = 0. 

zs difficult to deal with this equation because of the fractional Laplacian operator. But 
we can solve it via Poisson kernel. From [H], for x G (— r, r), 

sin g r (r 2 - x 2 )"/ 2 1 
18 



Figure 8. Escape probability in Example 14.21 b(x) = —x,a(x) = l,r = 2 



Obviously, p(—r) = 0. To justify p(r) = 1, we apply the substitution y = (r 2 — xv)(x — 
i;)" 1 to obtain 

sin zes r r 

p{r) = — / (r - u) a_1 (r 2 - v 2 )~%dv 



TT 

TTCt pi 



(1 — V) 2 V 2 1 d U 



Sm 2 
7T 

7T v 2' 2 ; 

^Xr ( f)r(i-f: 

71 2 I 

1, 



where the Beta and Gamma functions and their properties are used in the last two steps. 
The escape probability p(x) is plotted in Figured for various a values. 
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ff=0.25 




Figure 9. Escape probability in Example 14.31 r = 2 
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